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Abstract 

Nowadays in the MSSM the moderate values of tan /? are almost excluded by LEP II 
lower bound on the lightest Higgs boson mass. In the Next-to-Minimal Supersym- 
metric Standard Model the theoretical upper bound on it increases and reaches 
maximal value in the strong Yukawa coupling limit when all solutions of renormal- 
ization group equations are concentrated near the quasi-fixed point. For calculation 
of Higgs boson spectrum the perturbation theory method can be applied. We inves- 
tigate the particle spectrum in the framework of the modified NMSSM which leads 
to the self-consistent solution in the strong Yukawa coupling limit. This model al- 
lows one to get nih ~ 125 GeV at values of tan/? > 1.9. In the investigated model 
the lightest Higgs boson mass does not exceed 130.5 ± 3.5 GeV. The upper bound 
on the lightest CP-even Higgs boson mass in more complicated supersymmetric 
models is also discussed. 
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1 Introduction 



Last year there was a great progress in the Higgs boson searches. The experimental 
lower bound on the Higgs boson mass in the Standard Model (SM) has increased from 
95.2 GeV [0] to 113.3 GeV 0. At the same time the upper bound that comes from 
an analysis of radiative corrections to the electroweak observables has diminished to 
210 GeV P|. Thus the allowed region of the Higgs boson mass in the SM has shrunk 
drastically. Moreover at the LEPII a few additional bb events were observed 0]. They 
can be interpreted as a signal of the Higgs boson production with mass 115 GeV in the 
e'^e~ annihilation. Such Higgs boson mass does not agree with theoretical lower bound 
in the SM following from the stability of the physical vacuum up to the Planck scale 
Mpi ^ 2.4 X IQis GeV 

The simplest supersymmetric (SUSY) extension of SM is the Minimal Supersymmetric 
Standard Model (MSSM). Its Higgs sector includes two doublets Hi and H2. Each of 
them after spontaneous symmetry breaking acquires a nonzero vacuum expectation value 
Vi and V2, respectively. Instead of them, the sum of their squares v"^ = vf + f| and the 
value of tan/5 = V2/V1 are usually used. 

An important feature of supersymmetric models is the existence of a light Higgs boson in 
the CP-even Higgs sector. The upper bound on its mass strongly depends on the value 
of tan p. At the tree level the lightest Higgs boson mass does not exceed the Z-boson 
mass: < Mz \ cos2/?|. The loop corrections from the t-quark and its superpartners 
significantly raise the upper bound on m^: 



mh < yM|cos2 2/5 + A, 
where A in the one-loop approximation is given by 



_ 3 
^ 27r2 t;2 
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Here is the running top quark mass at the electroweak scale (at q = Mf°^^ = 174 GeV, 
Mf"^'^ is the t-quark pole mass), Xt is the stop mixing parameter, and Mg is the SUSY 
breaking scale which is expressed via the stop masses m^^ and m^^ • = y/^Wj^- The 
one-loop corrections @) attain maximal value at Xt = iv^M^. These corrections are 
proportional to and depend logarithmically on the SUSY breaking scale Ms- They 
are almost insensitive to the choice of tan/3. The absolute value of A is of order of M|. 
The one-loop and two-loop corrections to the lightest Higgs boson mass were calculated 
and analysed in and [jlO|, respectively. The upper bound on the hghtest Higgs boson 



mass grows with increasing of tan/? and ln(M|/m2), and for large tan/? (tan/? ^ 1) 
reaches 125 — 128 GeV. In the bounds on the mass of the Higgs boson in the SM and 
the MSSM were compared. 

However, the large values of tan (3 are unpreferable for two reasons. First of them is the 
proton decay. If one assumes that electroweak and strong interactions are embedded in 
the SUSY Grand Unified Theory (GUT) at high energies then too fast proton decay is 
induced due io d = b operators. The major decay mode is p — vK^ . The proton life-time 
Tp in this case is inversely proportional to tan^ /3. When tan (3 is large enough, calculated 
in the framework of SUSY GUT models proton life-time contradicts an experimental 
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restriction on it. Another problem of large tan/5 scenario concerns flavour-changing 
neutral currents. The branching ratio 6 — > 57 rises with tan/5 as Br{b —>■ 57) ~ tan^/3. 
Thus the values of tan /3 1 (tan /3 > 40) lead to the unacceptable large flavour-changing 
transitions. 

In the case of moderate values of tan/5 (tan/? < 5) the 6-quark and r-lepton Yukawa 
couplings are small and one can get an analytical solution of one-loop renormalization 
group equations [0 



For the top quark Yukawa coupling constant ht{t) and gauge 



couplings gi{t) the analytical solution has the following form: 

E{t) 

Yt{t) = 



QF{t) 



1 + 



QYt{0)F{t) 



ai{t) 



(3) 
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where Yt{t) 



ht{t) 



and aiit) 



l + 6i«,(0)t' 

t 

F{t) = j E{T)dT, 


, , , . The index i varies from 1 to 3, that 

An J ' ' \ Att J 

corresponds to U{1), SU{2), and SU{3) gauge couplings. The coefficients bi of one-loop 

beta functions of ai{t) are hi = 33/5, 62 = 1, &3 = —3. The initial conditions Yt{0) 

and ai{0) for the MSSM renormalization group equations are usually set at the Grand 

Unification scale Mx ~ 3 x 10^^ GeV where all gauge coupling constants coincide. The 

variable t is defined by a common way t = ln(M|./g2). 

Substituting the numerical values of the gauge couplings one finds that at the electroweak 
scale the second term in the denominator of the expression describing the evolution of 



Yt{t) is approximately equal to 



When /i^(0) > 1 the dependence of the top quark 



m?(o)' 

Yukawa coupling on its initial value Yt{0) disappears and all solutions are concentrated 
near the quasi-fixed point [^: 



^QFp(^o) 



(4) 



where to = 2\n{Mx/M^ 



polCN 



65. Together with Yt{t) the trilinear scalar coupling of 



the Higgs boson doublet H2 with stops At{t) and some combination of their masses 
9Jlj(t) = uiq + my + m\ are also driven to the infrared quasi-fixed points. In the vicinity 
of these points At{t) is proportional to the universal gaugino mass M1/2 at the scale 
Mx and OTt^ (t) ~ ^1/2- Although the solutions of the MSSM renormalization group 
equations achieve quasi-fixed points only for infinite values of l^t(O), the deviations from 

them at the electroweak scale are determined by the ratio — — — which is quite 



small if /i2(0) > 1. 



6F(to)rt(0) 



The behaviour of solutions of the MSSM renormalization group equations near the quasi- 
fixed point at tan /5 ~ 1 and particle spectrum have been studied by many authors [Il3|-|T5[] . 
It has been shown that in the vicinity of this point the h — r Yukawa coupling unification 



is realised [T^. In the recent publications (see pSHlH) the value of tan/? that corresponds 
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to quasi-fixed point regime has been calculated. It is restricted between 1.3 and 1.8. Such 
comparatively low values of tan (3 yield more stringent bound on the lightest Higgs boson 
mass in the MSSM. So, it does not exceed 94 ± 5 GeV [jl5|,|T6|. The obtained theoretical 
bound on nth has to be compared with the lower experimental one in the SM since it 
was computed for the SUSY breaking scale Ms of order of 1 TeV when all other Higgs 
bosons and superparticles are heavy enough. The straightforward comparison shows that 
the quasi-fixed point scenario and the considerable part of the MSSM parameter space 
are almost excluded by LEP II data that stimulates the theoretical analysis of the Higgs 
sector in nonminimal supersymmetric models. 

In this article the spectrum of the Higgs bosons in the Next~to-Minimal Supersymmetric 
Model is reviewed. The lightest Higgs boson mass in the NMSSM attains its maximum 
value in the strong Yukawa coupling limit, when the Yukawa couplings are much larger 
then the gauge ones. All supersymmetric models contain a large number of free parame- 
ters what is the main obstacle in the way of their investigations. For example, each SUSY 
model includes three or four independent SUSY breaking constants which determine the 
SUSY particles spectrum. Nevertheless, in the strong Yukawa coupling limit the solu- 
tions of the renormalization group equations are focused near the quasi-fixed point which 
simplifies the analysis. We propose a modification of the NMSSM which allows one to 
get rrih ~ 125 GeV for moderate values of tan/? and study Higgs boson spectrum of the 
model. In the last part the lightest Higgs boson mass in more complicated SUSY models 
is considered. 



2 Higgs sector of the NMSSM 

2.1 The /i— problem and parameters of the NMSSM 

The simplest extension of the MSSM is the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM). Historically the NMSSM was suggested as a solution of the /i-problem 



in the supergravity (SUGRA) models [Q. In addition to observable superfields these 
models contain a "hidden" sector where local supersymmetry is broken. In the super- 
string inspired SUGRA models the "hidden" sector always includes the singlet dilaton 
S and moduli Tm superfields. They appear in the four-dimensional theory as a result of 
compactification of extra dimensions. The full superpotential of SUGRA models can be 
presented as an expansion in powers of observable superfields 

W = Wo{S, T„) + T^){HiH2) + ht{S, T,){QH2)U'n + . . . , (5) 

where Wo{S,Tm) is the superpotential of the "hidden" sector. From the expansion it 
is obvious that parameter /i should be of order of the Planck scale because that is the 
only scale characterising the "hidden" (gravity) sector of the theory. On the other hand, 
if /i ~ Mpi then the Higgs doublets get huge positive masses m|^^ ~ /i^ ~ M|j and 
electroweak symmetry breaking does not occur at all. 

In the NMSSM a new singlet superfield Y is introduced. By definition the superpotential 
of this model is invariant with respect to the discrete transformations [19|. The 



symmetry usually arises in the superstring inspired models in which all observable 
superfields are massless in the exact supersymmetry limit. The term fi{HiH2) does 
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not satisfy the last requirement. Therefore it must be ehminated from the NMSSM 
superpotential. Instead of it the sum of two terms 

Wt, = \Y{hA) + (6) 



arises 



18|-P0| . After electroweak symmetry breaking the singlet field Y acquires a nonzero 
vacuum expectation value ((F) = y/\/2) and the effective /x-term (// = \y/\/2) is gen- 
erated. 

The NMSSM superpotential contains a lot of Yukawa couplings. But at the moderate 
values of tan/3 all of them are small and can be neglected except for /i^, A, and x. In 
addition to the Yukawa couplings the lagrangian of the NMSSM contains a large number 
of soft supersymmetry breaking parameters. Each of the scalar and gaugino fields has 
a soft mass (m^ and M,, respectively). Each of the Yukawa couplings corresponds to 
the trilinear scalar coupling Ai in the full lagrangian. The number of these unknown 
parameters can be considerably reduced if one assumes the universality of the soft SUSY 
breaking terms at the scale Mx- Thus, only three independent dimensional parameters 
are left: the universal gaugino mass M1/2, the universal scalar mass mo, and the univer- 
sal trilinear coupling of scalar fields A. Naturally universal soft SUSY breaking terms 



appear in the minimal supergravity model |21] and in the simplest models deduced from 



the superstring theories The universal parameters of soft supersymmetry breaking 
determined at the Grand Unification scale have to be considered as boundary conditions 
for the renormalization group equations that describe the evolution of all fundamental 
couplings up to electroweak scale or SUSY breaking scale. The complete system of the 
NMSSM renormalization group equations can be found in 



2.2 The CP— even Higgs boson spectrum 

The Higgs sector of the Next-to-Minimal Supersymmetric Standard Model includes six 
massive states. Three of them are CP-even fields, two are CP-odd fields, and one is a 
charged field. The determinants of the mass matrices of the CP-odd and charged Higgs 
bosons equal zero. It corresponds to the appearance of two Goldstone bosons: 

rf = ^2 sin /5 Im + ^2 cos /5 Im 
r/+ = sin 13 + cos P (H^)* 

which are swallowed up by the massive vector and Z-bosons during the spontaneous 
breaking of SU{2) ^U{1) symmetry. For this reason the masses of neutral CP-odd bosons 
and charged boson are easily calculated. 

In the CP-even Higgs sector the situation is more complex. The CP-even states arise 
as a result of mixing of the real parts of the neutral components of two Higgs doublets 
with the real part of the field Y. The determinant of their mass matrix does not vanish 
and in order to calculate its eigenvalues one has to diagonalize the (3 x 3) mass matrix. 
Instead of Re if}', Reifj; KeY it is much more convenient to consider their linear 
combinations: 

Xi = V2 cos P Re + V2smp Re H^, 

X2 = -V2 sin f3 Re H^ + V2 cos (3ReH^, (8) 
X3 = V2ReY. 
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In the basis (H) the mass matrix of CP-even Higgs fields can be simply written in the 



following symmetrical form (see [£5 



/ 9V 1 S'V d^V \ 



1 d^V 



V dvdf3 
1 d^V 



V dvd(3 f 2 d(3'^ 
d^V 1 d^V 



dvdy 
1 d^V 

V dydp 



(9) 



V dydj3 dy"^ j 



\ dvdy 

where V{vi,V2^y) is the effective potential of the NMSSM Higgs sector: 

V{Hi, H2, Y) = ml\Hi\^ + ml\H2\^ + ml\Y\^ + 
+ {\A^{H,H2)Y + ^A^Y^ + \x{H,H2){Y*f + h.c.) + 



(10) 



where /Sy{Hi, H2, Y) is the sum of loop corrections to the effective potential, g and g' 
are the gauge constants of the SU{2) and U{1) interactions, respectively {gi = ^/5/3g'). 

It is well known that the minimum eigenvalue of a matrix does not exceed its minimum 
diagonal element. Thus the lightest CP-even Higgs boson mass is always smaller than 



mi < M- 
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dv'^ 



v^ sin^ 2/3 + M| cos^ 2/5 + A , 
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In the right-hand side of inequality (|T3) A is the contribution of loop corrections to the 
Higgs boson potential. The expression (|ll|) was obtained in the tree level approximation 
(A = 0) in 1^ . The contribution of loop corrections to the upper bound on the lightest 
Higgs boson mass in the NMSSM is almost the same as in the minimal SUSY model. In 
particular, in order to calculate the corrections from the t-quark and its superpartners 
one has to replace the parameter /i in the corresponding formulas of the MSSM by Xy/^/2. 
The Higgs boson sector of the NMSSM and loop corrections to it were studied in p4|-p^. 
Let us also remark that the upper bound on rrih in the NMSSM was compared (see [0) 
with theoretical bounds in the SM and in its minimal supersymmetric extension. 

The calculation of CP-even Higgs spectrum is simplified in the most interesting realistic 
case when all superparticles are heavy {Ms ^ Mz)- In this case the contributions of new 



particles to the electroweak observables are suppressed as 



Ms 



(see, for example. 



On the other hand, the prediction for the values of strong coupling constant at the 
electroweak scale a^i^Mz) that can be obtained from the gauge coupling unification 
is improved with increasing of the supersymmetry breaking scale Ms- For Ms — 1 TeV 
it becomes close to the a^lMz) = 0.118(3) which has been found independently from the 
analysis of the experimental data [^]. Also it should be noted that the lightest Higgs 
boson mass reaches its maximal value in the SUSY models for Ms ~ 1 — 3 TeV. 



6 






0^ 






V12 




El 









V22 







El) 




V^31 


V32 


^33 



In the considered limit the mass matrix (P) has a hierarchical structure and can be 
represented as a sum of two matrices [^: 



Ml = I El I + I \/2i 1^22 1^23 I . (12) 

The first matrix is diagonal with E\ = ^ and El^ ^ Mg. The matrix elements Vu, V22, 
V33, V12 = V21 are of order of M|. The other matrix elements that correspond to mixing 
of Xi and X2 with xs are equal to 

^13 = ^31 = XvXu V2S = Vs2 = XvX2, (13) 

where Xi ~ X2 ~ Ms- 

Ml 

Considering the ratio as a small parameter it is easy to diagonalize the mass matrix 

-^2,3 

by means of usual quantum mechanical perturbation theory. For the Higgs boson 
masses it yields: 



-2 Y"^ 
2 ^ p2 I T/ I \2 2^'-i I \2 2 ^2 



77.2 ■ • ■ - T72 _ ■p2 
2 



4 - El + V22 + X'v^j;^. (14) 

E2 - E^ 

, ^ ^v^ sin^ 2(3 + Ml cos' 2(3 + A - X'^v'^. 
2 E^ 



2 



The explicit expressions for the Ef and Vij can be found in . For simplicity we restrict 
our consideration to the first order of perturbation theory and neglect matrix element 

M^ 

V\2 because its contribution to is of order of ^ 



Ml 



The perturbation theory becomes inapplicable when \El — E^\ \vX2. However, the 
mass matrix (12) can be easily diagonalised even in this case. In order to do this one 
should choose the basis where the matrix element M23 is zero. After that in the new 
basis the Higgs boson masses can be computed using ordinary perturbation theory. 



The first three terms in the last relation in ([1^ reproduce the upper bound on the lightest 
Higgs boson mass in the NMSSM. Their sum is equal to V\\ in our notations. The last 

Ml 



term in this expression gives a negative contribution to vtih. Even when the ratio 

goes to zero it does not vanish. Thus in the NMSSM the mass of the lightest CP-even 
Higgs boson can be considerably less than its upper bound . 



2.3 Renormalization of the Yukawa couplings and soft SUSY 
breaking terms 

According to inequality (PI) the upper bound on ra^ rises when A increases and the value 
of tan (5 diminishes. For tan /? 3> 1 the value of sin 2(5 goes to zero and the upper bound 
on the lightest Higgs boson mass in the NMSSM coincides with that one in the minimal 
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SUSY model. With decreasing of tan f3 the top quark Yukawa couphng at the electroweak 
scale /it (to) grows. The analysis of the solutions of the NMSSM renormalization group 
equations reveals that a rise of the values of the Yukawa couplings at the electroweak 
scale entails an increase of them at the Grand Unification scale. As a result the upper 
bound on the lightest Higgs boson mass in the NMSSM reaches its maximum value in 
the strong Yukawa coupling limit when the Yukawa couplings are much larger than the 
gauge ones at the scale Mx- 

The renormalization of the NMSSM coupling constants in the strong Yukawa coupling 



limit has been studied in p0|j3l| . In the considered case the Yukawa couplings are at- 
tracted towards a Hill type effective (quasi-fixed) line (x = 0) or surface (x ^ 0) which 
restrict the allowed regions of ht, A, and x. Outside this range the solutions of renormal- 
ization group equations blow up before the Grand Unification scale Mx and perturbation 
theory becomes invalid at ~ ^x- While the values of the Yukawa couplings at the 
scale Mx grow, the region, where all solutions are concentrated, shrinks and hf{0), A^(0), 
>i^(0) are focused near the quasi-fixed points These points appear as a result of 



intersection of the quasi-fixed line or surface with the invariant (fixed) line. The latter 
connects the stable fixed point in the strong Yukawa coupling regime with infrared 



fixed point of the NMSSM renormalization group equations ||3^. The properties of in- 
variant lines and surfaces were reviewed in detail in ||5|j3^. 



When the values of Yukawa couplings tend to quasi-fixed points the trilinear scalar 
couplings Ai[t) and some combinations of scalar particle masses Dn^{t), where 

ml{t)=ml{t)+ml{t)+ml{t), (15) 



become insensitive to their initial values A and StUq at the scale Mx ||3T| . For the universal 
boundary conditions one has 

A{t)=ei{t)A + f,{t)M^/2, 
w?^{t) = ai{t)ml + hit) Ml + Ci{t)AMy2 + di(t)A\ 

The functions ej(t), /j(t), ai(t), &j(t), Cj(t), and di{t) remain unknown since an analytical 
solution of the NMSSM renormalization group equations has not been found yet. While 
the Yukawa couphngs tend to infinity the values of functions ej(to), ^(^0); and diito) 
vanish. It means the solutions of renormalization group equations go to the quasi-fixed 
points too. In the vicinity of the quasi-fixed points Ai{t) are proportional to M1/2 and 

mm ~ Ml,,. 



3 Particle spectrum in the modified NMSSM 
3.1 The modified NMSSM 

The fundamental parameters of the NMSSM at the Grand Unification scale have to be 
adjusted so that the minimisation conditions of the effective Higgs boson potential are 
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satisfied: 



dV{vi,V2,y) ^ dV{vi,V2,y) ^ dV{vi,v2,y) ^ . . 

o „ U, U. l-l- ' / 

afi 0V2 oy 

Since the vacuum expectation value v is known they can be used for the calculation of 
A, mo, and Mi/2- But in the strong Yukawa coupling limit it is impossible to get the real 
solution of nonlinear algebraic equations (jl^)- Thus although the recent investigations 



35| , |36[| reveal that the upper bound on the lightest Higgs boson mass in the NMSSM 
is larger than the one in the MSSM by 7 — 10 GeV, in the considered region of the 
NMSSM parameter space the self-consistent solution cannot be obtained. Such solution 
of equations (|17|) appears for A^(0), x^(0) < 0.1 when the upper bound on rrih is the same 
as in the minimal supersymmetric model. 

Moreover due to the symmetry three degenerate vacuum configurations arise. After 
a phase transition at the electroweak scale the Universe is filled by three degenerate 
phases. The regions with different phases are separated from each other by domain walls. 
The cosmological observations do not confirm the existence of such domain walls. The 
domain structure of the vacuum is destroyed if the discrete symmetry of the NMSSM 
lagrangian disappears. An attempt of Z3 symmetry breaking by means of operators of 
dimension d = 5 was made in . It was shown that their introduction leads to quadratic 



divergences in the two-loop approximation, i.e., to the appearing of hierarchy problem. 
As a result the vacuum expectation value of Y turns out to be of the order of 10^^ GeV. 

Thus, in order to avoid the domain wall problem and get the self-consistent solution in 
the strong Yukawa coupling limit, one has to modify the NMSSM. The simplest way is 
to introduce the bilinear terms n{HiH2) and /i'F^ in the superpotential which are not 
forbidden by electroweak symmetry. At the same time one can omit the coupling x, that 
allows to simplify the analysis of the modified NMSSM. In this case (x = 0) the upper 
bound on the lightest Higgs boson mass reaches its maximum value. Neglecting all the 
Yukawa constants except for ht and A one gets the following expression for the modified 
NMSSM (MNSSM) superpotential ||]: 



l^MNSSM = l^{HiH2) + fi'V + XY{HiH2) + ht{H2Q)Ut (18) 



The bilinear terms in the superpotential (|T8]) break the Z3 symmetry and the domain 



walls do not arise, because the degenerated vacua do not exist. The introduction of 
the parameter permits to obtain the self-consistent solution of algebraic equations 
(p!7|) when /;,^(0), A^(0) ^ gfi^)- In the supergravity models the bilinear terms may be 
generated due to the additional terms [Z{HiH2) + Z'Y^ + h.c] in the Kahller potential 
3^,40] or due to nonrenormalizable interaction of the Higgs doublet superfields with the 



'hidden" sector ones 40.41 
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The effective Higgs boson potential of MNSSM can be written as tlie sum 



V{H,, H,, Y) = ^il\H,\' + fil\H,\^ + fil\Y\^ + 



+ XAxY{HiH2) + Xfi'Y*{H^H2) + A/iF \H^\^ + lifsP + h-c. 



+ 



V2 



+ \'\(H^H,)\' + X'Y^ \H,\' + \H,\']+^ \ - | + 



(19) 



where the parameters /z^ are expressed via the soft SUSY breaking terms as follows: 

2 2,2 2 2,2 2 2|/2 

fil = + fl , /i2 = m2 + /i , fly = TUy + fl , 

The /i-terms in the superpotential (^) lead to the appearance of bilinear scalar couplings 
B and B' in the effective potential of the Higgs fields. They arise as a result of the soft 
supersymmetry breaking. The values of B and B' depend on the mechanism of the fi 
and fi' generation. For a minimal choice of the fundamental parameters all soft scalar 
masses and all bilinear scalar couplings should be put equal at the scale Mx'- 

ml{Mx) = mliMx) = ml{Mx) = ml 
B{Mx) = B\Mx) = Bo. 



Thus in addition to the parameters of the SM the modified NMSSM contains seven 
independent ones: 

X,fi,fi',A,Bo,ml,Mi/2. 



3.2 The procedure of the analysis 

Although the parameter space of the considered model is enlarged it is possible to get 
some predictions for the Higgs boson spectrum in the strong Yukawa coupling limit. It 
is reasonable to start our analysis from the quasi-fixed point of MNSSM pO|: 



p?^P(to)= 0.803, pX'' {to) = 1-77, pg^^P(to) = 6.09, 
pfP(to)= 0.224, p5f (to) = -0.42, pg^f (to) = -2.28, 



(20) 



because all solutions of renormalization group equations are concentrated in the vicinity 
of it at hf{0), A^(0) ^ gfi^)- In the relations (pO|) the following notations are used: 

pt,x{t) = -^TTT, Yt = j-r^, Yx(t) = j—r, PA,. = and p^? - 



«3(t) ' ^ (47r)2' ^^ ^ (47r)2'^"'-^ My^ ''^ M 



2 • 

1/2 
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For given pt{to) (|20| ) the value of tan /5 can be extracted from the expression which relates 
the running t-quark mass mt{Mf°^^) with htito): 

mt{Mr'n = ^h,{Mr'> sin (3. (21) 

We substitute in the left side of equation (|2l|) the value of the running top quark mass 
mt(Mf°'*^) = 165 ± 5 GeV calculating in the MS scheme |^2[. The uncertainty in the 



mt{Mf°^'^) is determined by the experimental error with which the top quark pole mass 
is measured: Mf"^"^ = 174.3 ±5.1 GeV [^]. In the infrared quasi-fixed point regime that 
corresponds to maximal allowed values of /ij (0) and A^(0) (we take hf{0) = A^(0) = 10) 
we obtain tan/3 ^ 1.88 for mt{Mf"^'') = 165 GeV 



At the first step of our analysis the supersymmetry breaking scale is also fixed by means of 
the condition M3(1000 GeV) = 1000 GeV, where M3 is the gluino mass. This condition 
permits to calculate immediately the universal gaugino mass at the Grand Unification 
scale. It ensures that the superparticles are much heavier than the observable ones. 



Next we use the equations ([T7|) that define the minimum of the Higgs boson potential of 
the modified NMSSM to restrict the allowed region of the parameter space. Instead of fi 
it is more convenient to introduce /iefr = + ^y/V^- Then after some transformations 
we obtain: 



2 _ ml - ml tan^ /? + Az 1 , .2 
'•^s^" tan2/3-l 2^^^' 

XyX2 



sin 2/3 



-2 (^^+ cos 2/3; (22) 



ml + ml + 2filg + fv^ + 
X 

Y, _ A 



y (ml + /.'^ + 5V) = ^v'X^ - A. 



where Aj are the contributions of loop corrections and 

Xi = (2/iefr + (/i' + Ax) sin 2/3) , X2 = (/i' + Ax) cos 2/3. 

As the values of v and tan (3 are known one can find from the equations (^) the vacuum 
expectation value y and parameters Bq and /iefr- In the numerical analysis we take into 
account only the loop corrections from t-quark and its superpartners, because they give 
a leading contribution. Therefore Aj are the functions of and do not depend on Bq 
and y. From the first equation of (|2^ ) the absolute value of fies is calculated. The sign of 
Heff is not determined and should be considered as a free parameter. The bilinear scalar 
coupling Bq and vacuum expectation value y are computed from the two other equations 
of (|2^). The last of them points out that the value of y is of the order of Xv'^/Ms and 
much smaller than v if the superparticles are heavy enough. 

Since /iefr, Bq, and y have been found we investigate the dependence of the Higgs boson 
spectrum on A, mg, and /i' using the relations ([TBI) . For the masses of CP-odd states 
one can get an exact analytical result: 
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m\ = ml + ml + + y t;^ + Aa, 

2 2 I /2 D' ' I 2 , A 

mQ = my + ii - B ^ + —V +A3. 



(23) 



The mass matrix of CP-even Higgs sector has a hierarchical structure and can be written 



in the form (|T2|) with 



El 
^11 



m\ + ml + 2/i^g, 



M| cos^ 2/3 + -A^t;^ sin^ 2/3 + An, 



V^12 = V, 

^33 
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lAV-^M|)sin4/3 + Ai2, 



M| sin^ 2/3 + -^v' cos^ 2/3 + A^ + A22, 



V31 = At;Xi + Ai3, 
1 
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V32 = At;X2 + A 



23; 



-X'v' + A33. 



(24) 



In the formulas (^) and ( ^4]) Aq, A3, A^ and Ajj (An = A) are the loop corrections 
to the CP-odd and CP-even mass matrices. The mass matrix of CP-even Higgs sector 
can be diagonalised and the expressions for the masses of CP-odd states ( p3D can be 
simplified using the perturbation theory of quantum mechanics. In the main order of 
perturbation theory the masses of the heavy Higgs bosons are: 

m^, and m^^ 
similar to given in 



El-, 



ml 



El 



while the first order perturbation corrections have the form 



3.3 Numerical results 

The results of the numerical analysis of the particle spectrum near the MNSSM quasi- 
fixed point are presented in Figs. 1-3. There are two regions of the MNSSM parameter 
space. In one of them the mass of the lightest CP-even Higgs boson is larger than the 
upper bound on ruh in the MSSM (see Figs, la and 2a) whereas in the other region it is 
smaller (see Figs, lb and 2b). As follows from the relations (|l^ the lightest Higgs boson 
mass in the NMSSM and in its modification attains its upper bound when V13 (or Xi) 
goes to zero. In the MNSSM it happens if 

, _ 2/ieff V2A13 
^ ""iha5"^"" At;sin2/3- ^^^^ 

Thus mh is larger in that part of the parameter space where the signs of /x and /x' are 
opposite. If /i' tends to infinity then the singlet CP-odd and CP-even fields get huge 
masses and their contribution to the effective potential of the Higgs bosons vanishes due 
to the decoupling property. In the considered limit the lightest Higgs boson mass is the 
same as in the minimal SUSY model. 

It is necessary to emphasize that the one-loop and even two-loop corrections give an 
appreciable contribution to the mass of the lightest CP-even Higgs boson. So the two- 
loop corrections |l^ reduce its mass approximately by 10 GeV. They nearly compensate 
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the growth of hghtest Higgs boson mass with increasing of SUSY breaking scale Ms 
which arises because of one-loop corrections. Due to loop corrections the values of rrih 
for /ioff > and /iefr < become different. The contribution of loop corrections A rises 
as the stop mixing parameter Xt = At + /Xefr/ tan/3 increases. Since near the quasi-fixed 
point v4j < (see ([20|)), the absolute value of Xt is larger if /ieg < 0. 



As the part of the parameter space where /ieg and /i' have the same signs is almost 
excluded by Higgs searches data from LEP II, we investigate the particle spectrum in the 
case when the signs of /Xgff and /x' are opposite. In the most interesting region, where 
the hghtest Higgs boson mass is close to its upper bound, the value of /i' is considerably 
larger than /Xeff and Ms (see (^)). Moreover the product B' n' is positive. Indeed from 



the second relation of the system (22), which defines the minimum of the MNSSM Higgs 
boson potential ([T9|) , it follows that the bilinear scalar coupling B and have different 
signs. As a consequence near the maximum of the curves in Figs, la and 2a the sign of 
B' coincides with the sign of /i'. 

It means that the heaviest particle in the modified NMSSM is the CP-even Higgs boson 
that corresponds to the neutral field Y . Its mass is m| > /i'^ and is substantially larger 
than the scale of supersymmetry breaking. As one can see from Figs. Ic and 2c the mass 
of the other heavy CP-even Higgs boson itih is almost insensitive to the value of /z' since 
m| 3> m?^. The masses of CP-odd states are always smaller than /x'^. If the value of ^' 
grows the mass of the heaviest CP-odd state increases too as ~ /i'^. When the 
value of yu' diminishes the lightest CP-odd boson mass m^a decreases and for /i' ~ B' 
becomes of the order of electroweak scale. At given values of /i' the low constraint on 
yu' appears which comes from the requirement of m^^ > 0. However even if the mass of 
the lightest CP-odd state is of the order of M^, it will be quite difficult to observe it in 
future experiments because the main contribution to its wave function gives the CP-odd 
component of the singlet field Y . The heaviest fermion in the modified NMSSM is the 
neutralino {frty^^) which is a superpartner of the scalar field Y . Its mass is proportional to 
/i'. The remaining masses of neutralinos (^«yJ, charginos (m-±), squarks, and sleptons 
do not depend on the /i'. 

The spectrum of new fermion states, squarks and sleptons is also insensitive to the choice 
of the parameter A, since near the quasi-fixed point ( pO|) the dependence of scalar masses 
ml and trilinear scalar couplings Ai on it disappears. For this reason the lightest CP- 
even Higgs boson mass is almost independent of A. Nevertheless the dependence of the 
heavy Higgs boson spectrum on the universal trilinear scalar coupling A is conserved. It 
occurs because the bilinear scalar coupling B' is proportional to A. The dependence of 
the Higgs boson masses on the parameter A for mo = is presented in Figs. Id and 2d. 

Although everywhere in Figs. 1 and 2 we put mo = 0, the qualitative pattern of the 
particle spectrum does not change if the universal soft scalar mass varies from zero to 
Mfj^. It should be noted that the masses of squarks, sleptons, heavy Higgs bosons, heavy 
charginos and neutralinos rise with increasing of mo while the spectrum of the lightest 
particles remains unchanged. 

Up to now the particle spectrum in the quasi-fixed point regime which corresponds to 
the initial values of the Yukawa couplings /i^(0) = A^(0) = 10 has been studied. In 
the vicinity of the quasi- fixed point (H) for mt(MP°''') = 165 GeV and M3 < 2 TeV 
the lightest Higgs boson mass does not exceed 127 GeV. The results presented in Table 
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point out that the quahtative pattern or MNSSM particle spectrum does not change 
even if h^{0) ^ A^(0) or h'^{0) <^ A^(0) as long as the Yukawa couplings are much larger 
than the gauge ones at the scale Mx- The masses of the Higgs bosons and the masses 
of the superpartners of the observable particles were calculated there for the values of 
/i' computed by the formula (|25| ) and A = mo = 0. At the same time from the Table 
one can see that the numerical value of the lightest Higgs boson mass is raised from 
105 - 113 GeV for X\0) = 2 to 118 - 128 GeV for X'^{0) = 10. 

Therefore at the last stage of our analysis we investigate the dependence of the upper 
bound on nih on the choice of the Yukawa couplings at the Grand Unification scale. 
For each h^(tQ) we find the value of tan/3 using the relation (^I]) and choose A^(to) and 
/i' so that the lightest Higgs boson mass attains its upper bound. The obtained curve 
m/i(tan /?) is plotted in Fig. 3, where the upper bound on mh in the MSSM as a function 
of tan/9 is also presented. Two bounds are very close for large tan/? (tan/? ^ 1). The 
curve m/j(tan/?) in the MNSSM reaches its maximum when tan/? = 2.2 — 2.4 which 
corresponds to the strong Yukawa coupling limit. The numerical analysis reveals that 
the mass of the lightest CP-even Higgs boson in the modified NMSSM is always smaller 
than 130.5 ± 3.5 GeV, where the uncertainty is mainly due to the error in the top quark 
mass. 



4 Final remarks and conclusions 

We have argued that the upper bound on the lightest Higgs boson mass in the NMSSM 
attains its maximal value in the strong Yukawa coupling limit. In the considered limit all 
solutions of renormalization group equations gathered near the quasi-fixed points. If the 
scale of supersymmetry breaking is much larger than the electroweak one the perturbation 
theory can be used for the calculation of the Higgs boson masses. However even when 
Ms ^ Mz the lightest CP-even Higgs boson mass in the NMSSM is appreciably smaller 
than its upper bound in the dominant part of the parameter space. Besides in the strong 
Yukawa coupling limit within the NMSSM with a minimal set of fundamental parameters 
the self-consistent solution does not exist. Moreover the Z3 symmetry of the NMSSM 
superpotential leads to the domain wall problem. 

We have suggested such a modification of the NMSSM that allows to get the self- 
consistent solution in the strong Yukawa coupling limit and at the same time to avoid 
the domain wall problem. The superpotential of the modified NMSSM (MNSSM) in- 
cludes the bilinear terms which break the Z3 symmetry. We have studied the spectrum 
of the Higgs bosons in the MNSSM. The qualitative pattern of the particle spectrum is 
most sensitive to the choice of two parameters - /i' and Ms- The limit n' ^ Ms, when 
the CP-even and CP-odd scalar fields become very heavy, corresponds to the minimal 
SUSY model. In the most interesting region of the parameter space, where the mass of 
the lightest Higgs boson is larger than that in the MSSM, the Higgs boson mass matrix 
has the hierarchical structure and can be diagonalised using the method of perturbation 
theory. The heaviest particle in this region of the MNSSM parameters is the CP-even 
Higgs boson corresponding to the neutral scalar field Y and the heaviest fermion is Y 
which is the superpartner of the singlet field Y. The lightest Higgs boson mass in the 
considered model may reach 127 GeV even for the comparatively low value of tan /? ~ 1.9 
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and does not exceed 130.5 ± 3.5 GeV. 



The obtained upper bound on the mass of the hghtest CP-even Higgs boson is not 
an absolute one in the supersymmetric models. For instance, the upper bound on rrih 
increases if new 5 + 5 supermultiplets appear in the NMSSM at the SUSY breaking 
scale. These multiplets change the evolution of gauge couplings. Their values at the 
intermediate scale rise if a number of new supermultiplets increases. For this reason the 
allowed region of the Yukawa couplings at the electroweak scale is expanded. It leads to 
the growth of the upper bound on ruh as the number of 5 + 5 supermultiplets rises (see 
Fig. 4). The investigations performed in |3^ showed that the introduction of four or five 



5 + 5 supermultiplets raises the theoretical bound on the lightest Higgs boson mass up to 
155 GeV. If more than five multiplets are introduced at the SUSY breaking scale then 
the gauge couplings blow up before the Grand Unification scale and perturbation theory 
is not valid at ~ ^x- 

Recently the upper bound on the lightest Higgs boson mass in more complicated SUSY 
models has been analysed P^-^ . In particular, in addition to the gauge singlet superfield 



three SU (2) triplets Tj with different hypercharges can be introduced into the Higgs boson 
superpotential: 

W^Higgs = XYiHA) + XiiHifoH^) + XiiHifA) + X2(^2Tii^2) + • ■ • • (26) 
As a result the expression for the upper bound on rrih changes 



ml < Ml cos^ 2(3 + 



+ M ) sin2 2/5 + 2x1 cos^ /? + 2x1 sin' P 



+ A. (27) 



The appearance of triplet superfields destroys the gauge coupling constant unification at 
high energies. In order to restore the unification scheme of the electroweak and strong 
interactions one has to add several S'[/(3) multiplets, for example four (3 + 3), which do 
not participate in the SU{2) ® U{1) interactions. A numerical analysis |45|j46[] reveals 
that the unification of gauge couplings then occurs at the scale Mx ~ 10^'' GeV. As one 
can see from Fig. 5 (see also |^) the upper bound on the lightest Higgs boson mass rises 
with growth of tan/5 and for tan/5 ^ 1 can be approximately equal to 190 GeV |H6|. 



Also the upper bound on mh is raised if in the MSSM a fourth generation of the quarks 



and leptons exists ^^l- However up to now there has not been found any evidence of 
the existence of the fourth generation in the SM or the MSSM. Moreover, new particles 
give considerable contributions to the electroweak observables which upset the agreement 
between theoretical predictions and the results of experimental measurements. Thus the 
growth of the upper bound on the lightest Higgs boson mass in the supersymmetric 
models is usually accompanied by substantial increase in the number of particles that 
may be considered as the main drawback of these models. 
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Table. Spectrum of the superparticles and Higgs bosons for mt(174 GeV) = 165 GeV, 
A — mo — and for different initial values of /i^ (0) and A^(0) (all mass parameters are 
given in GeVs). 
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Figure captions 

Figure 1. The particle spectrum (in GeVs) in the modified NMSSM as a function 
oi z = /i'/lTeV and x = A/My2 for h'^iQ) = X^{0) = 10, ml = 0, M3 = ITeV, and 
/ieff < 0. Thick and thin curves in Fig. la and lb correspond to the hghtest Higgs boson 
mass calculated in the one- and two-loop approximation, respectively. In Fig.lc and 
Id thick and thin curves reproduce the dependence of CP-even Higgs boson masses ms 
and rriH on z and x while dotted and dashed curves represent the CP-odd Higgs boson 
masses rriAi and as a functions of these parameters. The dashed-dotted curve in 
Fig.lc corresponds to the mass of the heaviest neutralino. 

Figure 2. The particle spectrum (in GeVs) in the modified NMSSM as a function 
of z = /i'/lTeV and x = A/M1/2 for /^^(O) = X^{0) = 10, = 0, M3 = ITeV, and 
/ieff > 0. The notations are the same as in Fig.l. 

Figure 3. Upper bound on the lightest Higgs boson mass (in GeVs) in the MSSM 
(lower curve) and in the modified NMSSM (upper curve) as a function of tan/? for 
Ms = 2 TeV. 

Figure 4. The dependence of the upper bound on the lightest Higgs boson mass (in 
GeVs) in the NMSSM (solid and long-dashed curves) and in the MSSM (short-dashed 
curve) on the value of tan j3. The solid curve is the upper bound on rrih in the NMSSM 
with four additional 5 + 5 multiplets while the long-dashed curve is that one in the 
ordinary NMSSM (see ||). 

Figure 5. Theoretical upper bound on the lightest Higgs boson mass in the supersym- 
metric model, which contains singlet and three SU{2) triplet fields in the Higgs sector 
and four 3 + 3 colour multiplets, as a function of tan/? for Ai = (see [0). 
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